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Intensive larval samplings of mosquitoes (Diptera: Culicidae) were carried out 
between 2005 and 2008 in several biotopes located in a varying climate region in 
Eastern Spain. The biodiversity was analyzed and divided into alpha, beta and 
gamma components with the aim of comparing the mosquito biodiversity ac- 
cording to the different structure of the landscape due to the incidence of climatic 
and anthropic patterns. Likewise the synanthropic index of Nuorteva was calcu- 
lated for each species. A total of 11,279 mosquitoes belonging to 29 species was 
collected and identified. Mosquito biodiversity is higher in the wettest and non- 
anthropized areas. Using a cluster analysis, all this information was also used to 
group the different regions studied depending on its mosquito fauna. Moreover 
the re-emergence of antroponosis, like malaria, seems unlikely given the low val- 
ues of the synanthropic index for the anophelines captured. 


R. Bueno Mari & R. Jiménez-Peydro, Entomology and Pest Control Laboratory, 
Cavanilles Institute for Biodiversity and Evolutionary Biology, University of Va- 
lencia, Official Post 22085, 46071 Valencia, Spain; E-mail address of the corre- 
sponding author: ruben.bueno@uv.es 


Received 29 July 2010, accepted 10 January 2011 


1. Introduction 


Mosquitoes, like other organisms show a spatial 
variation in relation to several factors such as en- 
vironmental heterogeneity or habitat and host 
preferences (Zhong et al. 2003). Understanding 
the links between habitats, environmental factors 
and occurrence of immature mosquitoes (Dipte- 
ra: Culicidae) is essential for an efficient applica- 
tion of mosquito control methods (Pemola Devi 
& Jauhari 2007). 

To compare mosquito diversity according to 
the structure of the landscape, the division of the 
alpha, beta and gamma components is considered 


to be a very useful tool (Whittaker 1972). The al- 
pha diversity (œ) is the specific richness of a com- 
munity that we consider homogeneous. The beta 
diversity (D) refers to the degree of replacement in 
the specific composition between different com- 
munities of a landscape. And lastly we can define 
the gamma diversity (y) as the specific richness of 
the grouped communities that form a landscape, 
resulting from both alpha and beta diversities. 
This method of biodiversity analysis is useful not 
only to explore the climatic, physical or biologi- 
cal influences on biodiversity, but also to study 
the effects of human pressure on biodiversity 
(Halffter 1998, Moreno 2001). Although it is a 


ENTOMOL. FENNICA Vol. 22 * Mosquito biodiversity in Spain 191 


method much used to estimate the biodiversity of 
many insects, there is little information about its 
use in mosquito studies. 

Our aim was to investigate the changes in 
mosquito biodiversity due to environmental (cli- 
matic and anthropic) factors in Eastern Spain. 


2. Material and methods 
2.1. Study area 


In Spain studies of mosquito biodiversity are 
scant and outdated, mostly concerned with times 
when malaria was endemic (Bueno Marí & Jimé- 
nez Peydró 2008). A total of 64 mosquito species 
have been reported in Spain, but some of them, 
like Anopheles labranchiae Falleroni or Aedes 
aegypti (Linnaeus), are currently considered 
eradicated or absent (Bueno Marí & Jiménez 
Peydró 2011). Consequently, records of only 54 
species are considered valid at present (Eritja ef 
al. 2000). For our study we selected the Valencian 
Autonomous Region (Fig. 1) due to its heteroge- 
neity of climate. In this Mediterranean area we 
differentiated between five regions or sectors ac- 
cording to their patterns of temperature and pre- 
cipitation (GVA, 2003): 


— Setabense sector (S): Characterized by a typi- 
cal Mediterranean climate, with low oscilla- 
tions of mean temperature (« 13°C) between 
winter and summer, and peak precipitations in 
autumn. 

— Valenciano-Tarraconense sector (VT): Simi- 
larly represents a typical Mediterranean cli- 
mate, but with higher homogeneity in the spa- 
tial distribution of precipitations than S. 

— Maestracense sector (MS): Characterized by 
a Mediterranean climate strongly influenced 
by a Continental climate. There are wide tem- 
perature oscillations (« 18*C) between winter 
and summer, and peak precipitations during 
spring. Moreover it is the coldest and wettest 
region of the study area. 

— Manchego sector (MN): Despite it also has a 
Continental influence, it can be easily distin- 
guished from MS due to its lower topography. 

— Alicantino-Murciano sector (AM): Charac- 
terized by a dry Mediterranean climate, show- 





Jf 


Fig. 1. Situation of the study area with the five coro- 
logical sectors. AM: Alicantino-Murciano sector, MN: 
Manchego sector, MS: Maestracense sector, S: Seta- 
bense sector and VT: Valenciano-Tarraconense sec- 
tor. 


ing an intermediate climate between typical 
Mediterranean and Desertic climate. 


To study the effects of human pressure on mos- 
quito biodiversity, each larval biotope was typi- 
fied according to its anthropization degree as fol- 
lows: 


— Urban biotope or high anthropization area: 
situated in permanent human settlements or 
surroundings. 
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Table 1. Indices and formulae used to estimate alpha, beta and gamma diversities. 








Estimation Indices Formulae 

Alpha diversity Margalef (Dug) Dy, = S-1t/in N 
Simpson (A) A= p? 
Shannon-Wiener (H) H'--Xp,ln p, 


Equity of Pielou (J) 


Jaccard (I) 
Whittaker (£) / (B)* 


Gamma (y) 


Beta diversity 


Gamma diversity 


J’=H’/ Hmax-H'IIn S 


l'-clatb-c 
B,- Sla-1l(By-1-d 


y"axfxe 





# S: species richness, N: total number of individuals, p; frequency of species i, a: number of species present in site A, b: number 
of species present in site B, c: number of species present in both site A and B, d: average number of samples occupied by one 
species, a: average number of species in a community, e: total number of communities studied. 


* Modification proposed by Schluter & Ricklefs (1993) to calculate the average f. 


— Ruralbiotope or medium anthropization area: 
placed in temporary human settlements (live- 
stock or agricultural areas) or surroundings. 

— Wild biotope or low anthropization area: po- 
sitioned in natural areas without apparent hu- 
man actions around it. 


This habitat characterization allowed us to apply 
the Nuorteva synanthropic index (Nuorteva 
1963) to evaluate the domiciliation degree of 
each species caught: 
(S) - Qa;* b, —2c)/2 (1), 
where a, b and c represents the percentage of cap- 
tures ofthe i species in urban, rural and wild areas 
respectively. The values of the index range from 
+100 (highest degree of domiciliation) to —100 
(lowest degree of domiciliation). The domicili- 
ation degree refers to the preference of different 
species to breed and/or develop in anthropic envi- 
ronments. From the point of view of disease vec- 
tors, the domiciliation degree is postulated as a 
very important evolutionary issue to explain the 
emergence and endemism of human diseases that 
previously were exclusive zoonoses. The Nuor- 
teva index was initially proposed for the study of 
synanthropic flies, but has also been used with 
culicid mosquitoes (Forattini et al. 1993). Given 
the fact that human settlements are frequently 
linked to the worsening of water quality, the use 
ofthis index could be interesting, not only to eval- 
uate the anthropophilic degree of each species, 
but also to elucidate their tolerance levels of water 
eutrophication. 


2.2. Sampling methods 


Using simple random sampling, we sampled lar- 
val forms from multiple sites within each study 
sector using the dipping method (Service 1993) 
for 8 months a year (March-October) 2005—2008 
(4 years). Data were collected from all identifi- 
able aquatic environments across 23,260 km of 
the study area. 


2.3. Data analysis 


Mosquitoes were identified according to the keys 
of Encinas Grandes (1982), Darsie & Saminadou 
Voyadjoglou (1997) and Schaffner et al. (2001). 
The nomenclatural criteria were based on Reinert 
(2000). Data were analyzed statistically using the 
Biodiversity.R package (Kindt & Coe 2005). For 
each environment we calculated several biodiver- 
sity indices (Table 1). To calculate the ecological 
distance between different environments, we de- 
veloped a cluster analysis and a principal compo- 
nent analysis (PCA), both based on Jaccard dis- 
tance. The cophenetic correlation was also calcu- 
lated for each Jaccard cluster with the aim of cal- 
culating the degree of reliability ofthe classifica- 
tion system used (R Development Core Team 
2005). 


3. Results 


In the larval captures of 679 different sampling 
points, a total of 11,279 mosquitoes belonging to 
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Table 2. — A. Number specimens of species in five genuses captured in different environments, and synanthropic 
indices (SI). — B. Estimators of alpha biodiversity. For the abbreviations of corological sectors, see Fig. 1. 











Corological sectors Areas 

AM MN MS S VT Urban Rural Wild SI 
A. No. of specimens and SI 
Aedes 
A. vexans 0 89 0 8 85 12 6.5 
A. vittatus 0 0 0 1 50 0 5 46 -79 
Anopheles 
A. algeriensis 16 9 0 20 0 0 33 12 12.5 
A. atroparvus 0 4 61 0 26 0 25 66 —56.5 
A. claviger 0 7 25 2 0 0 2 32 -79 
A. maculipennis 0 0 61 0 13 0 9 65 -67 
A. marteri 0 0 3 0 0 0 0 3 —50 
A. petragnani 0 7 147 187 145 0 88 398 —61 
A. plumbeus 0 0 171 0 0 0 171 0 50 
Culiseta 
C. annulata 0 0 4 2 35 0 18 23 -37 
C. longiareolata 239 102 693 788 797 411 1,379 829 13 
C. subochrea 1 2 10 8 108 60 69 -35.5 
Culex 
C. hortensis 18 44 455 75 111 0 207 496 —58.5 
C. impudicus 0 7 30 22 26 0 20 65 —64 
C. laticinctus 1 47 205 177 482 158 451 303 —12.5 
C. mimeticus 7 15 74 123 123 0 111 231 —53.5 
C. modestus 2 0 5 122 70 0 164 35 17 
C. pipiens 587 61 346 1,395 1,129 729 2,011 778 21 
C. territans 10 27 197 133 121 0 158 330 —51 
C. theileri 26 2 0 6 42 0 50 26 0.5 
Ochlerotatus 
O. berlandi 0 0 42 0 0 0 42 0 50 
O. caspius 312 12 4 103 207 69 364 205 21 
O. detritus 60 0 0 24 77 13 90 58 2 
O. echinus 0 0 45 89 11 0 119 26 0.5 
O. geniculatus 0 0 4 49 9 0 48 14 6.5 
O. gilcolladoi 0 0 2 13 0 0 0 15 —50 
O. pulcritarsis 0 0 14 0 0 0 14 0 50 
Orthopodomyia 
O. pulchripalpis 0 0 10 0 0 0 10 0 50 
Uranotaenia 
U. unguiculata 0 0 0 0 28 0 17 11 0.5 
B. Alpha diversity 
Abundance 1,279 346 2,087 3,339 3,618 1,380 5,751 4,148 
Specific richness (S) 12 14 24 20 21 5 27 25 
Margalef index (D,,,) 1.54 2.22 2.91 2.34 2.44 0.55 3 2.88 
Simpson index (A) 0.69 0.84 0.87 0.76 0.83 0.62 0.8 0.88 
Shannon-Wiener index (H) 1.44 2.09 2.37 1.9 2.19 1.14 2.17 2.46 
Equity of Pielou index (J) 0.58 0.79 0.75 0.63 0.72 0.36 0.66 0.89 





29 species were identified (Table 2). The compar- 
ison of a biodiversity indices reveals that the MS 
sector is the most diverse (Dye = 2.91, à = 0.87, 
H’ = 2.37) and the AM the least diverse Ou. = 
1.54,2 =0.69, F’ = 1.44). MS sector is the wettest 


sector; it features bodies of water in greater quan- 
tity and typological diversity. 

Concerning the degree of environmental 
anthropization, we suggest that human settle- 
ments negatively influence mosquito biodiversi- 
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Table 3. Estimators of beta biodiversity for each community comparison. For the abbreviations of corological sec- 
tors, see Fig. 1. The average £ for comparisons of corological sectors is 0.32 and that for areas is 0.51. 





Comparisons Jaccard (/)) 


Jaccard distance (Cluster) 


Whittaker (B, ) 





Corological sectors 
AM-MN 

AM-MS 

AM-S 0.6 
AM-VT 0.5 
MN-MS 
MN-S 
MN-VT 
MS-S 
MS-VT 
S-VT 


Areas 

Urban-Wild 0.2 
Urban-Rural 
Wild-Rural 0.8 


"— —— —— — — 
O Q io Q 0 t I 0 T w 
IDOCOZONAND 


o 
N 
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o © 





ty: urban biotopes (e = 0.55, A = 0.62, H' = 
1.14), rural biotopes (D4,73.4 —0.80, H’=2.17) 
and wild biotopes (Dy = 2.88, A = 0.88, H’ = 
2.46). With the exception of tree hole breeding 
mosquitoes, associated mostly with forest distri- 
bution, the synanthropic index of Nuorteva indi- 
cates that anophelines, excepting Anopheles 
algeriensis Theobald (S = 12.5), show low 
domiciliation degree, whereas Culex pipiens 
Linnaeus, Ochlerotatus caspius (Pallas), Culiseta 
longiareolata (Macquart) and Culex modestus 
Ficalbi are species better associated with habitats 
which present surrounding human settlements. 
Furthermore, other species, like Culex hortensis 
Ficalbi (S — —58.5), Culex impudicus Ficalbi (S = 
—64), Culex territans Walker (S — —55) and Culex 
mimeticus Noè (S ^ —53.5), exhibit a strong pref- 
erence for wild environments. 

The analysis of f biodiversity (Table 3) indi- 
cates that S and VT are the closest sectors in their 
specific compositions (J, = 0.71), the specific re- 
placement among them being consequently very 
low (B. = 1.07). It is important to remember that 
both sectors are the only representatives of the 
typical Mediterranean climate in our study area. 
With regard to the degree of environmental 
anthropization, the rural and wild biotopes are the 
closest (7, — 0.80), and consequently also show a 
low specific replacement (B. = 1.16). 

The dendrograms elaborated from Jaccard 


AM 


MS 
0.5 


0.3 


0.70 Urban 


0.65 


0.60 


0.55 


Wild Rural 
Fig. 2. Cluster analysis using Jaccard distance. — a. 
Corological sectors, cophenetic correlation (r,) = 0.98. 
— b. Degree of anthropization, (r.) = 0.88. For the 
abbreviations of corological sectors, see Fig. 1. 
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Fig. 3. Principal component analysis based on Jaccard distance. For the abbreviations of 
corological sectors, see Fig. 1. PC values (PC1, PC2) of species within the boxes: Upper box: 
Aedes vittatus (0.2, 0.2), Anopheles algeriensis (—2.9, 0.5), A. atroparvus (-3.2, —3.1), A. clavi- 
ger (-2.7,-0.9), A. maculipennis (—3.1, —3.1), A. marteri (-0.2, —0.2), Culiseta annulata (—0.1, 
—0.2), C. subochrea (—0.2, —0.1), Culex impudicus (-0.2, —3.1), C. mimeticus (0.2, —2.8), C. mo- 
destus (0.5, 0.4), C. theileri (—0.1, —0.2), Ochlerotatus berlandi (—2.9, —2.6), O. echinus (-0.2, 
—1.1), O. geniculatus (—0.2, —0.8), O. gilcolladoi (—0.1, —0.7), O. pulcritarsis (—0.2, —0.7), Ortho- 
podomyia pulchripalpis (—0.2, —0.6) and Uranotaenia unguiculata (—0.2, —0.5); Lower box: Culex 
territans (0, —8) and Anopheles plumbeus (-3.3, —8.1). 


distance (Fig. 2) corroborate the same conclu- 
sions already given. The high values of jaccard 
distance cophenitic correlations (imas = 0.98, 
ae 0.88) indicate high correlation between 
the ecological distance observed in our study and 
the distance predicted in the hierarchical configu- 
ration of the clusters. 

The PCA analysis shows that C. hortensis and 
O. caspius are the species that best define the axes 
closest to MS (wettest and highest biodiversity) 
and AM (driest and lowest biodiversity), respec- 
tively (Fig. 3). The opposite nature of the axes 
seems to indicate that the species should also ex- 
hibit opposite behaviours. 

In accordance with the values of the 
synanthropic index, the second PCA also reveals 
that C. hortensis, C. mimeticus, C. territans and 


Anopheles petragnani Del Vecchio best define 
the axis closest to wild environments (Fig. 4). 

Finally, y (climatic) 18 29.12 and y ontopic) 18 29.07. 
These both are virtually identical to the value of 
the total species richness (29 species) found in the 


study area. 


4. Discussion 


Results from our study show that mosquito diver- 
sity varies across climates and land use patterns. 
Previous researchers have proposed the use of 
mosquitoes as bio-indicators of forest degrada- 
tion (Dorvillé 1996) and of anthropic pressure on 
other natural environments (Montes 2005). MS is 
the wettest sector, consequently the abundance of 
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Fig. 4. Principal components analysis based on Jaccard distance. PC values (PC1, PC2) of 
species within the box: Aedes vittatus (0.8, —1), A. vexans (1.6, —2.1), Anopheles algeriensis 
(71.1, 0.2), A. atroparvus (—2.8, —3.1), A. claviger (2.2, —0.8), A. maculipennis (1.9, —1), A. mar- 
teri (—0.4, 0.2), Culiseta annulata (—2.2, —1.7), C. subochrea (—2, —1.5), Culex impudicus (0.3, — 
0.5), C. modestus (1.5, —2.3), C. theileri (-0.8, —1.1), Ochlerotatus detritus (2.2, —2.1), O. berl- 
andi (0.1, 0.2), O. echinus (—0.2, 0.6), O. geniculatus (0.3, 0.9), O. gilcolladoi (—0.5, —0.9), O. 
pulcritarsis (0, —0.4), Orthopodomyia pulchripalpis (0.4, —0.6), Uranotaenia unguiculata (1.6, —2). 


water and hetereogeneity of aquatic environ- 
ments increase the appearance of colonizable 
habitats for many species, thereby increasing 
mosquito biodiversity. On the other hand, AM is 
the driest sector, thus only a few species were ca- 
pable of adapting to this adverse situation. Simp- 
son indices show dominance by one or more spe- 
cies that are characterized by their high propor- 
tional abundance. In general, regardless of the 
sector, C. pipiens and C. longiareolata are the 
dominant species. In addition, high dominance 
should also be noted in C. hortensis in MS, O. 
caspius in AM and Culex laticinctus Noé in VT. 
Our study reveals that the wettest regions present 
higher mosquito biodiversity than the driest ones, 
especially in Mediterranean climates where tem- 
perature is not a limiting factor. 

The values of Nuorteva's synanthropic index 
are well related with the host preferences and tol- 
eration of different water qualities described for 


the different species catched (Schaffner et al. 
2001). However, in this situation there is proba- 
bly also an underlying human factor, since the 
coastal areas are the regions where human settle- 
ments tend to accumulate most, while inland 
areas have higher mosquito biodiversity. This hy- 
pothesis is also supported by the data related to 
the biodiversity in urban, rural and wild environ- 
ments. In this respect it is well known that C. 
hortensis is a zoophile (mainly feeding on ba- 
trachians and reptiles) and orophilic species, 
which prefers to colonize freshwater environ- 
ments and can complete its development in very 
varied breeding sites (Schaffner et al. 2001, 
Bueno Mari et al. 2009a). In contrast, O. caspius 
is a species characterized by a high anthropo- 
philic behavior, very common in coastal areas of 
the Paleartic Region and highly associated with 
temporary puddles of brackish waters (Rioux 
1958, Schaffner et al. 2001). In conclusion, a high 
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anthropization does not imply a reduction of mos- 
quito population in urban environments, but 
rather a selection filter that is only surpassed by a 
few species. 

Moreover, Fig. 4 indicates that the co-exis- 
tence of C. hortensis, C. mimeticus, C. territans 
or A. petragnani could be used as a bioindicator 
element for wild environments. On the other 
hand, C. pipiens and C. longiareolata are species 
that have a high biological plasticity and there- 
fore can breed in environments with widely vary- 
ing degrees of human pressure. 

Faunistically, it is important to note that ofthe 
29 mosquito species observed in this work, 13 
have never been recorded previously in the study 
area (Romeo Viamonte 1950, Encinas Grandes 
1982). This study allowed the capture of species 
poorly and irregularly found in Spain such as 
Anopheles marteri Prunelle and Uranotenia un- 
guiculata Edwards (Torres Cafiamares 1979, 
Bueno Mari et al. 2010), among others. The 
synanthropic indices for the anophelines confirm 
the low malariogenic potential described for 
Spain (Bueno Mari & Jiménez Peydró 20102). 
Moreover, no A. labranchiae were found, the 
Spanish distribution of which has previously 
been circumscribed south of our study area. The 
dramatic reduction of irrigated agriculture in 
Southeastern Spain perhaps led to the eradication 
of A. labranchiae (Eritja et al. 2000). This species 
was an important malaria vector and its capture 
was common up to 1946 (Clavero & Romeo 
Viamonte 1948). Since then it has not been re- 
corded in Spain (Blázquez & De Zulueta 1980, 
Bueno Marí & Jiménez Peydró 2010b). 

Of concern is the expanding distribution of 
Aedes albopictus (Skuse), the Asian tiger mos- 
quito. This species, potential vector of several 
arboviruses such as Dengue, Yellow Fever or 
Chikungunya among others (Bueno Marí & Ji- 
ménez Peydró 2010c), was first recorded in Spain 
in 2004 (Aranda et al. 2006) and was recently de- 
tected south of our study area (Bueno Marí et al. 
2009b). 
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